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Circulation control is a viable active flow control approach that can be used to meet the NASA 
Subsonic Fixed Wing project’s Cruise Efficient Short Take Off and Landing goals. Currently, 
circulation control systems are primarily designed using empirical methods. However, large 
uncertainty in our ability to predict circulation control performance has led to the development of 
advanced CFD methods. This paper provides an overview of a systematic approach to developing 
CFD tools for basic and advanced circulation control applications. This four-step approach includes 
“Unit”, “Benchmark”, “Subsystem”, and “Complete System” experiments. The paper emphasizes 
the ongoing and planned 2-D and 3-D physics orientated experiments with corresponding CFD 
efforts. Sample data are used to highlight the challenges involved in conducting circulation control 
computations and experiments. 


Nomenclature 


AOA 

= 

angle of attack (degrees) 

w = 

slot width (inches) 

AR 

= 

aspect ratio 

X = 

streamwise tunnel coordinate (inches) 

b 

= 

span (inches) 

Y = 

horizontal cross plane tunnel coordinate (inches) 

CC 

= 

circulation control 

Z = 

vertical cross plane tunnel coordinate (inches) 


= 

pressure coefficient 



Cd, 

= 

drag coefficient 

Symbols: 

C d 

= 

section drag coefficient 

a = 

pitch angle of attack (deg) 

C L 

= 

lift coefficient 

5 = 

jet angle relative to free stream (deg) 

Q 

= 

section lift coefficient 

5* = 

boundary layer displacement thickness 

Cm 

= 

momentum coefficient 

y = 

ratio of specific heats 

C, c 

= 

chord, wing, airfoil (inches) 

p = 

density 

h 

= 

slot height (inches) 

4> = 

separation angle relative to model (deg) 

LE 

= 

leading-edge 

r = 

circulation 

TE 

= 

trailing-edge 

e = 

separation angle relative to free stream (deg) 

MAC 

= 

mean aerodynamic chord (inches) 



M 

= 

Mach number 

Subscripts 

m 

= 

mass flow rate (lbm/sec) 

J = 

jet 

NPR 

= 

nozzle pressure ratio 

o = 

free stream total condition 

P 

= 

pressure (psf) 

00 — 

free stream static condition 

q 

= 

dynamic pressure (psf) 



r 

= 

radius of Coanda surface or flap (inches) 


Re 

= 

Reynolds number 



R 

= 

gas constant 



S 

= 

planform area (in 2 ) 



T 

= 

temperature (°F) 



U 

= 

velocity component in streamwise direction (ft/sec) 

V 

= 

velocity component in horizontal cross plane (ft/sec) 

w 

= 

velocity component in vertical cross plane (ft/sec) 
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Introduction 

Recent interest in circulation control (CC) technology has been prompted by the National Aeronautic 
and Space Administration (NASA) Cruise Efficient Short Take Off and Landing (CESTOL) initiative and 
the Air Force Research Laboratory (AFRL) Advanced Joint Air Combat System program. These initiatives 
highlight the need for improved high-lift or powered-lift systems to meet takeoff and landing goals that are 
integrated into transonic cruise configurations. Because circulation control technology will require 
significant modifications to the aircraft design it will need to “buy” its way onto the aircraft. This will 
require designers to optimize the entire wing and propulsion system with particular attention on air delivery 
system efficiency, safety concerns related to an engine out condition, and control issues related to large 
aerodynamic moments. 

The rapid turnaround required by designers to complete trade studies normally limits the use of 
Computation Fluid Dynamic (CFD) tools to the improvement of empirical methods. Empirical methods 
derived from experiments are traditionally used in the design process. However, the complexities of many 
flow control techniques in use today, including circulation control, far exceed the experimental databases 
used to develop the empirical methods. By coupling physics-based CFD tools with validation experiments 
we can develop parametric studies that characterize the sensitivities of wing geometries and blowing 
systems, required to advance empirical techniques necessary to design a CESTOL type aircraft. 

The application of the state of the art CFD tools to a CC airfoil with an experimental database was 
described at the 2004 NASA/ONR Circulation Control workshop. 1 Results from that workshop showed 
that CFD simulations were unable to consistently 
and accurately predict the performance of a CC 
airfoil. An example of one of the CFD predictions 
is shown in Figure 1. The failure to match the jet 
separation and the streamline turning for this 
example was linked to turbulence models 11 ’ 111 and 
CFD grid issues. In addition to identifying 
inconsistencies in CFD prediction capability, the 
workshop also identified that the experimental 
databases were inadequate for CFD validation. 

This paper describes a systematic approach to 
developing experimental and computational 
databases for improving CC prediction capability. In 
general, CFD validation is defined by determining 
how well the CFD model predicts the performance 
and flow physics when used for its intended 
purposes. lv The level of CFD validation can be 
defined by the complexity of the code and the experiment being used for validation, as described in Figure 
2. v These levels of validation are being pursued using NASA’s airfoil, semi-span wing, and full-span 
Hybrid Wing Body geometries that are the focus of this paper. 

For the purposes of this paper the term circulation control implies using a tangential jet on a highly 
curved aerodynamic or Coanda surface. The aerodynamic characteristics of circulation control wings have 
been experimentally and numerically studied since Henri Coanda’ s near fatal crash in 1910 vl ’ vll,vm Many of 
these studies have concentrated on the lift performance related to trailing edge shape, slot height, and 
blowing rates. Despite the large potential performance improvements related to circulation control, only 
two fixed wing flight demonstrations lx,x,xl have been evaluated. The results of both demonstrator programs 
showed large benefits of circulation control applications for short runways. Without requirements for short 
runway access, however, circulation control applications were not pursued and the remaining work was 
largely a laboratory effort. 

The renewed interest in CESTOL aircraft at NASA has been prompted by the escalating airport 
congestion and noise abatement criteria. xll,xl11 The FAA X1V has projected that air traffic will increase by a 
factor of 2 or more within the next 20 years. Airport noise has been identified by the FAA as the largest 
problem for 50 of the nation’s airports. xv One way to alleviate these problems is to use the shorter runways 
and the infrastructure that are already in existence at under utilized small airports. xvl,xvn It may also be 
necessary to alter the traffic patterns and glide slopes to keep the noise confined to the airport 
boundaries. xvlll,xlx These issues have prompted engineers at NASA to evaluate the improvements in aircraft 
take off and landing capabilities that can be integrated into efficient cruise geometries. 



Figure 1. Comparison of experimental and 
Reynolds Average Navier Stokes (RANS) CFD 
lift performance data for a typical CC airfoil 
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NASA has set goals xx,xxl for a state-of-the-art 100 passenger airliner to be operational by the year 2022 
and to have the following characteristics: 

• take-off and landing distance of less than 2,000 feet, 

• cruise Mach number >0.8, 

• 1,400 to 2,000 mile range capability, 

• noise containment within an airport footprint, and 

• low speed maneuverability 

Improvement and optimization of advanced flow control systems and wing geometries are required to 
achieve these goals. Several advanced flow control and propulsion technologies are being considered to 
meet these goals, including CC, upper surface blowing, over the wing vectored blowing and/or 
combinations of these technologies. To achieve an optimized CESTOL geometry, it is necessary to 
evaluate the requirements of the aircraft. xxll,xxm These current CESTOL goals are being continually 
influenced and refined by issues such as fuel prices and new air traffic requirements. However, pursuing 
these goals using CC may provide a low speed performance capability that exceeds what can be achieved 
with conventional high lift systems. In the remainder of the text we will discuss the use and validation of 
CFD for advanced CESTOL type of aircraft. 

Code Validation Process for Circulation Control 

Circulation control applications are difficult to compute reliably using state-of-the-art CFD methods as 
demonstrated by the inconsistencies in CFD prediction capability described in the 2004 NASA/ONR 
Circulation Control workshop. 1 For Reynolds-averaged Navier-Stokes methods, some turbulence models 
such as Menter SST k-omega xxlv and Spalart-Allmaras with rotation and curvature correction (SARC) XXV 
have been shown to predict CC flows reasonably well only at certain conditions while other methods have 
been shown to be very sensitive to numerical parameters. There is also a tendency for the predictive 
capability of CFD to degrade as the blowing increases. While it is possible that turbulence models bear 


Unit 

Problems 

* Simple Geometry 
Hardware 

* One Element of 
Complex Flow Physics 

* One Relevant Flow 
Feature 


* Low Experimental 
Uncertainty 

* All Initial Conditions 
and Boundary Conditions 
Measured 


Benchmark 

Cases 

• Special Hardware 
Fabricated 

* Two Elements of 
Complex Flow Physics 
■ Two Relevant Flow 
Features 


• Moderate Experimental 
Uncertainty 

■ Most Initial Conditions 
and Boundary Conditions 
Measured 


Subsystem 

Cases 

■ Subsystem or 
Component Hardware 

■ Moderately Complex 
Flow Physics 

* Multiple Relevant Flow 
Features 


* Large Experimental 
Uncertainty 

■ Some initial Conditions 
and Boundary Conditions 
Measured 


Complete 

System 

■ Actual System 
Hardware 

■ Complete Complex 
Fiow Physics 

* Multiple Relevant Flow 
Features 


* Limited Experimental 
Data 

* Most Initial Conditions 
and Boundary Conditions 
Unknown 



K7" 


3 



Jet Exit 

Velocity Profile 


2 -D CC Airfoil w/ Upper and 
Lower Blowing 



3-D CC Semi-span 
Wing w/ LE & TE 
Blowing 



3-D Hybrid Wing Body w/ 
USB augmented w/ CC Witig 


Increasing Geometric And Flow Complexity 


Figure 2. Four levels of CFD validation used to study circulation control 
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some responsibility, some of the disagreement may be due to an increasing loss of two-dimensionality in 
nominally 2-D experiments as the blowing rates increase. Unfortunately, it is difficult to determine the 
sources of the discrepancies because most of the CC experiments have not been adequately documented for 
CFD validation. 

One of the parameters that CC 
performance is typically characterized by and 
that must be carefully documented is the 
thrust or momentum coefficient, Cp. The 
sensitivity of the airfoil performance to Cp is 
dependent on the jet characteristics and the 
airfoil geometry,™ 1 particularly on the 
surface near the jet exit. There are two 
physical regimes that define circulation 
control as a function of blowing.™ 11 These 
regimes are commonly referred to as 
separation control and super-circulation 
control and exhibit different global 
efficiencies as determined by the change in 
unit lift due to the change in unit blowing as 
shown by the data in Figure 3. Figure 3. Circulation control flow regimes showing 

The physical description of the sensitivities and influence of wall blowing (WB) 

efficiencies of these regimes is demonstrated by the relationship of the jet separation location on the trailing 
edge surface and the interaction of the jet with the on-coming flow. Nominally the flow separation on the 
flap or Coanda surface is a function of the pressure gradient along the surface. As the thin jet is applied, the 
flow separation is moved aft along the surface. This jet also entrains and accelerates the flow in the 
vicinity of the wall-bounded jet, thus increasing the turning of the streamlines near the trailing edge. This 
also affects the leading edge stagnation region thus increasing the overall circulation around the airfoil and 
the lift. 

The transition from one regime to another is not always clearly identified and is dependent on the 
sharpness of the trailing edge. The systematic study of this mechanism should include different trailing 
edge geometries that include a hinged flapped geometry (i.e. fixed separation located at tip of flap) and a 
circular geometry (i.e. separation location free to move up to 180°) as shown in Figure 4. 

NASA’s approach to resolve the 
issues of CFD inconsistency and the 
inability to match experiments at 
high-blowing conditions involves a 
strong interaction between CFD and 
experimental investigations. For 
CFD it is believed that most of the 
challenges exist in establishing 
consistent requirements for 
boundary conditions, along with 
guidelines for numerical 
methodology, such as the use of 
sufficient grid refinement as well as 
transition and other parameters that can influence the solution. The planned development of CFD tools for 
circulation control applications utilizes RANS and LES codes. The benefits of understanding the physics 
of complex flow fields to improve the prediction and modeling techniques in CFD codes should also be 
used to improve appropriate experimental methods. 

The goal of this series of investigations is to acquire error bounded experimental results suitable for 
CFD validation with no attempt to achieve maximum performance on a specific geometry. The approach 
for these experiments is to develop a 2-D fundamental circulation control database that will be followed by 
a 3-D database that captures realistic physics to advance predictive tool development. 



(a) Flap Configuration (b) Circular configuration 

Figure 4. Example of two Coanda surfaces. 
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An often-overlooked aspect of the quality of 
information from CC experiments pertains to the 
quantification of the blowing parameters. An example 
of this is quickly shown by considering the methods 
used to determine Cjli. The momentum coefficient can 
be determined by the quantities shown in either 
Equation 1 or 2. Equation 1 expresses the momentum 
coefficient in terms of basic geometric and flow 
parameters. For small-scale experiments, the 
uncertainty of the slot height (h) dominates Equation 1 , 
because the absolute slot height measurement error (e.g. 

±0.002”) is typically large relative to the total slot 
height. This measurement can also be complicated by 
non-uniformities in the slot along the airfoil span and by 
small changes in slot height when the internal duct is 
pressurized. Using Equation 2 simplifies the required 
measurements for determining Cpby eliminating the 
need to measure slot geometry or jet density. These 
parameters are integrated into the mass flow rate, which 
can be measured directly. An additional common error is in the calculation of the jet velocity, Uj. A 
commonly accepted method for calculating the jet velocity assumes that the jet flow expands isentropically 
to the freestream static pressure as shown in Equation 3. This introduces an error since the static pressure 
at the jet exit is less than the freestream static pressure (as demonstrated in surface pressure measurements) 
at the jet exit. When Equation 1 is used, the error is exacerbated since the Uj value is squared. While the 
momentum uncertainties are typically smaller when using Equation 2, it is often a CFD validation 
requirement that the individual parameters be measured. The uncertainties for these parameters are 
cumulative and are dependent on the transducers and the instrumentation used for the measurement.™ 111 ’™ 
It is recommended that multiple measurement techniques be used to quantify these uncertainties. 


THRUST 2hb JET p JET Uj ET 

q<x>S Cb Po0 Ij2 


C P = 


THRUST mUj ET 


qooS qooS 

Where, 


( 2 ) 


Viet = 

(3) 


2'i'bjET(T 0 j E T ) 


.7-1 


y — 1 


P 

"|_ 00 


P 

V r °JET 2 


and 


m = PjETUjEl(hX W ) 


(4) 


Description of Validation Cases 

In this section, we will provide a description of the experiments in terms of the systematic approach 
shown in Figure 2. Example experimental and CFD results will emphasize the critical features related to 
circulation control physics. 


Unit Problem 

One of the main factors influencing the CFD solution is the 
characterization of the wall jet at the jet exit. This can be 
set with a boundary condition of specified velocity or 
momentum either at the jet exit plane, or at an inner 
plenum wall. In the latter case, the flow (including 
turbulence) is allowed to adjust naturally as it is forced out 
of the high-pressure plenum. An example showing the 
complexity of the results from a computation that includes 
a plenum is shown in Figure 5. In the former case, it is 
more difficult to prescribe realistic profiles at the exit 
plane. Unfortunately, both methods have traditionally 
relied upon under-resolved details from the experiment. 
The only information available from experiments 
regarding the jet has often been the jet momentum 
coefficient, Cp, as determined by Equation (1). In this 
equation, the jet velocity is typically obtained from 
conditions inside the plenum, combined with isentropic 
flow relations. This methodology introduces additional 
uncertainty into the CFD simulation. Measurements of 
flow conditions, especially the velocity at the jet exit, help 
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Figure 5. Flow near the upper-surface 
blowing slot of an airfoil, showing 
normalized turbulent viscosity contours 
along with mean flow velocity vectors (the 
plenum is to the left of the lower channel) 
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to remove this uncertainty by establishing definitive boundary conditions to match at and near the jet exit. 

Implementation details of the CFD boundary condition at the plenum inflow are probably not too 
important, as long as the plenum is of sufficient size and provides the desired momentum flux. For many 
studies, the velocity and density (i.e., momentum) are specified at the upstream face inside the plenum, and 
the pressure is extrapolated from the interior of the domain. Other internal plenum boundary conditions 
were also tried by the authors, such as specifying total pressure along with total temperature (and 
extrapolated pressure), but it made little difference to the solution. 

The small size of the jet slot height often 
complicates measurement of the jet profile. For 
small scale wind tunnel experiments, the slot 
height can range from 0.010” < h < 0.60”. A 
nominal slot height of 0.020” is typical of 
geometries described in this paper. Miniature 
pitot probes have an outer diameter of 0.010” 
resulting in a large bias due to integration in the 
large velocity gradients typical of wall bounded 
jets. The use of hot wire probes does not come 
without measurement difficulties that include 
interpretation of transonic hot wire results. The 
hot wire prongs can be as large as 5% to 10% of 
the slot height as shown in figure 6. This is 
often complicated by model and probe 
vibration. All of these factors lead to 
measurement uncertainties that must be 
quantified. 

Another relevant factor related to the jet exit flow is the oncoming boundary layer above the lip of 
the slot. Because this boundary layer interacts/mixes with the jet, it is important that it be characterized 
correctly in the CFD solution. Accurate experimental measurements of boundary layer profiles and 
turbulence properties in this region are needed for validation. 

Benchmark Cases 

Two benchmark cases will be described as part of the systematic CC validation effort. Two of the 
desired features for each of the benchmark CFD validation studies are a simple geometry that is capable of 
super-circulation and measured critical boundary conditions required by modern CFD codes. Even though 
the geometries are simple, a 2-D lift coefficient of 7 - 8 can be achieved in the super-circulation regime. 

The 2-D Benchmark case 

Two-dimensional CFD has proven itself to be useful for many aerospace design and trade studies. 
However, when validating turbulence models for specific new applications, it is extremely important to 
ensure that consistent conditions (geometry and boundary) are being considered for the computation and 
the experiment. For CC airfoils, this can be difficult to ascertain because of the sensitivity of the Coanda 
surface flow separation to jet boundary conditions and other factors that can be problematic to measure, as 
well as the inherent difficulty of maintaining two-dimensionality in the experiment at high blowing 
conditions. In ongoing work at NASA we are attempting to address these issues, with the ultimate goal of 
validating the capabilities of existing turbulence models for 2-D CC flows over a wide range of blowing 
conditions. 

To reduce uncertainty and to optimize the measurement capability, two independent wind-tunnel tests 
were conducted using the same model. The NASA LaRC Basic Aerodynamic Research Tunnel (BART) 
and the Georgia Tech Research Institute’s Model Test Facility (MTF) are comparable in size and speed, but 
emphasize different measurement techniques. The BART test series emphasizes the external flow physics 
by identifying the leading edge stagnation point, determining the jet separation location, and determining 
the jet trajectory for selected blowing conditions. The MTF test series characterizes the CC model’s high- 
lift blowing performance as a function of slot height, angle of attack, and Cp. 



Velocity Magnitude, ft/s 


Figure 6. CFD velocity profile at the jet exit with 
scaled hot wire probe, Cp=0.119, NPR=1.2, h=0.020” 
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PICOLO TUBES 


ELLIPTIC FOREBODY 


TE SLOT HEIGHT 
0.005" <h <0.045” 



LE SLOT HEIGHT* 
0.010” <h <0.025” 
(PHASE 2) 


Figure 7. NASA/GTRI 
CC-E0020EJ airfoil 


r/C=9.44% 


2-D CFD Validation Model, 


A 2-D geometry with an elliptic 
leading edge and a large circular trailing 
edge profile, shown in Figure 7, was 
selected for the 2-D validation 
experiment. This geometry was chosen as 
the benchmark geometry for NASA’s 2-D 
effort because it is characteristically 
simple and has a large trailing-edge radius 
for accurate measurements of jet 
separation. The dual blowing (both upper 
and lower surfaces) capability will be 
used to manage the drag characteristics, 

but will have a lower priority in the initial validation database experiments. 

While it has been difficult to compute 2-D circulation control applications reliably, it has also been 
difficult to generate 2-D experimental data independent of the influences from wind tunnel walls. The two- 
dimensionality of the flow and the wall effects are proportional to the amount of circulation control or lift 
generated by the wind tunnel model. The comparison made in Figure 1 between the two-dimensional CFD 
simulations and the wind tunnel measurements show how CFD tends to over predict the performance of a 
circulation control model at high lift coefficients where tunnel wall effects become large. These effects are 
generally related to the walls limiting the streamline turning, as shown in Figure 8, and the wall and model 
juncture for 2-D testing shown in Figure 9. 


TUNNEL BOUNDARY 



Figure 8. 2-D CFD prediction of the floor and 
ceiling walls, AR=3.26, Cp=0.115, NPR=1.2 Q=4.87 
(no walls, blue), Q=4.43 (walls present, red), 
h/C=0.0023, wind tunnel height to chord=4.88, CC- 
E0020EJ airfoil 



WIND TUNNEL 
WALL JUNCTURE 


Figure 9. 3-D CFD prediction of the wall 
juncture flow, AR=3.26, Cp=0.23, NPR=1.4 
Q=5.09, h/C=0.0023, CC-E0020EJ airfoil 


The effect of the wind tunnel walls in the 2-D plane can be estimated by comparing the performance 
results with and without floor and ceiling walls. The presence of the walls decreased the lift coefficient by 
approximately 10% compared to free air calculations for the case shown in Figure 8. Three dimensional 
effects due to interactions with the spanwise walls are difficult to characterize in 2-D experiments. In 
similar experiments e.g., Novak et al.,™ the interaction of the sidewall boundary layers with the wing and 
its Coanda jet produced vortical structures that induced downwash along the span of the wing, yielding 
relatively large recommended angle-of-attack corrections (e.g., Aa = -8.94 deg for Cp = 0.226). Estimated 
corrections such as this are necessary when comparing to 2-D CFD computations. If sidewall blowing or 
other sidewall boundary layer control is used, it is helpful to quantify its influence on the flowfield at 
various blowing rates. 
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An example of the Particle Image Velocimetry (PIV) data compared to 2-D RANS calculations, shown 
in Figure 10, illustrates how CFD can be used to estimate the AOA effects of the wind tunnel walls. The 
leading-edge stagnation location measured with the PIV data is at x/c value of 4.03%. The comparable 
CFD calculation in free air (i.e. no wind tunnel walls) identifies the stagnation to be at an x/c of 7.75%. An 
AOA adjustment of -5° for the CFD simulation resulted in a close match to the experimental results for 
both the stagnation location and airfoil pressure distribution. The AOA adjustment based on classic wall 
corrections to account for tunnel wall interference is of this magnitude. Further investigation of wall 
interference using CFD tools is warranted. 




Figure 10. Example of PIV stagnation streamline data compared to CFD computations, GACC 
airfoil at (Xgeometric = 0° 


The measurements 
of the external flow field 
are driven by two CFD 
validation requirements, 

(1) capture the global 
streamline behavior 
related to blowing, and 

(2) measure the boundary 

layer and jet 
characteristics on the 
airfoil. The 

measurement suite for 
this series of tests was 
driven by model scale 
and optical access. 
External flow 

measurements will be 
obtained using a 2-D PIV 



system. This system will capture the details of the external flow to include the edge of the boundary layer. 
PIV measurements near the surface are problematic at this scale. Boundary layer details have been 
measured using a single-sensor hot wire. Examples of PIV measurements for a similar experiment are 
shown in Figure 11. 

The boundary layer transition characteristics for low-Reynolds number testing are often difficult to 
assess. For this series of experiments, the tunnel velocity was increased until constant baseline drag was 
achieved. Based on this criteria, a minimal tunnel operational limit was established for the validation 
experiments of q=15 psf, corresponding to Re c =0.50xl0 6 . Due to the presence of a laminar separation 
bubble near the leading edge of the model, a transition strip was added to the lower surface at x/c=3%. 
Example pressure distributions are shown in Figure 12. 

Because accurate computations of the jet are crucial, it is particularly important to measure the as-built 
slot width as a function of span to ensure 
adherence to specifications and minimal 
spanwise variation. Based on experience, the 
degree to which the slot geometry changes 
under pressure should be determined. In terms 
of boundary conditions, the jet conditions need 
to be measured as well as the freestream 
conditions at some specified location upstream 
in the wind tunnel. 

In terms of the desired measurements for 
CFD validation, surface pressures (including 
spanwise variations) are necessary, along with 
integrated force and moment data. Boundary- 
lay er/jet-^ velocity and turbulent shear stress 
profiles are the most useful quantities for 
judging the capabilities of the turbulence 
models. Without these measurements, it is 
difficult to make assessments because 
sometimes CFD values agree reasonably well 
with experimental surface pressures for the wrong reasons. In addition, the measured boundary-layer 
profiles on the airfoil upstream of the jet slot are extremely important to make sure that the “external” flow 
interacting with the jet is the correct one. Accurate measurements of the jet separation location off the 
Coanda surface as well as the stagnation streamline location near the leading edge are also key measures 
used to qualify the success of a CFD simulation. Imaging in the far-field away from walls can be useful, 
but it is not as important for the CFD validation as the other quantities listed above. 

Characterization of the boundary layers (including jet) is generally needed at Reynolds numbers 
typical of most of the wind tunnel tests. In other words, where do the boundary layers and jet become 
turbulent? Is there a laminar bubble present on the airfoil? Is the jet issuing from the slot laminar, 
turbulent, or transitional? 

Semi-span Benchmark Experiments 

The primary objective of this experiment is to evaluate Reynolds number effects for advanced 3-D 
circulation control wings for both high-lift and cruise configurations using the same model. This unique 
data set will also provide a benchmark database for CFD validation. Many of the issues described for the 
2-D benchmark validation experiment are still applicable for the semi-span experiment. These benchmark 
tests are being planned for the NASA FaRC National Transonic Facility (NTF) and the NASA FaRC 14- 
by-22-Foot Subsonic Tunnel to begin in the summer of 2009. 

The design criteria for this series of experiments is based on a realistic cruise configuration that 
contains the blowing capability for circulation control powered lift and cruise enhancement. To be 
consistent with NASA CESTOF system studies™ a swept wing with an aspect ratio of 9.5 is desired. 
However, achieving this aspect ratio for the NTF semi-span study would require scaling the slot to 
unreasonably small values. The compromise was to design the outboard section of the wing for NTF 
(AR~5) and extend the wing to hybrid wing configuration for future low-speed testing in the FaRC 14-by 
22-Foot Subsonic Tunnel and high-speed testing at NASA Ames Research Center’s 11-Foot Pressure 
Tunnel. 



X/C 

Figure 12. Experimental pressure distributions and 
corresponding lift for the CC-E0020EJ 
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While the goals of this semi-span CFD validation study are not intended to optimize wing geometry, it 
is desired to choose a state-of-the-art geometry with characteristics that are consistent with the CESTOL 
goals. Selection of the wing profile was therefore driven by the cruise performance (Mach > 0.8), thus 
limiting the trailing edge geometry to internally blown flaps. 

The NASA TMA-0712 supercritical profile selected for the 3-D validation experiments will be 
modified to accommodate an upper- surface blowing slot that will be located at the break line of the trailing 
edge hinged flap shown in Figure 13. This wing will be fabricated using a modular approach. The model 
will have a removable trailing edge hinged flap and a removable leading edge that can accommodate 
leading edge blowing or variations with and without slats, or droop, etc. It is intended that the wing can be 
used for other approaches to active flow control for future development as well. 

The NTF test will take advantage of cryogenic and pressure capabilities to achieve high Reynolds 
number conditions. Projected tunnel conditions will be limited to -50°F and 5 atmospheres, resulting in a 
Re c = 15xl0 6 at Mach=0.2 and Re c =30xl0 6 at Mach=0.8. The air delivery system for NTF is being 
developed for advanced blowing concepts that include 3-D circulation control. Volumetric flows must be 
controlled for a wide range of Cp’s, as shown in a typical operational map in Figure 14. A Cp or NPR 
sweep would typically operate at an average slot height ,e.g., h=0.060 average shown with red line in 
Figure 14. The maximum jet flow condition will be limited to a jet mach number of 1. The air delivery 
and mass flow systems must be capable of low flow and high flow operations. These conditions typically 
require a dual flow control system that includes temperature control. 


Two Independent Air Supplies 
(Directed to TE orLE or Combination) 



Figure 13. Sketch of NTF circulation control semi-span TMA-0712 model 


The conceptual design of the semi-span circulation control model has dual flow paths for independent 
control of the two air delivery systems. Internal pressure vessels are being designed to be consistent with 
the circulation control operational requirements described above. Model design and fabrication are being 
performed in parallel with the development of the new NTF air delivery system. Initial testing is planned 
for the 3rd quarter of 2009. 

The CFD modeling challenges discussed in the 2-D CFD development are noticeably increased for 
even a simplistic three-dimensional 
configuration such as that proposed for the 
NTF experiments. The increased flow 
field complexity is not the only pacing 
item, the growth in the grid complexity and 
density required to properly model the 
internal flow passages is another 
significant factor. Due to wind tunnel 
model size constraints and typically high 
plenum operating pressures, internal 
support structures are often required to 
maintain model integrity and geometric 
accuracy. Thus, a realistic CFD simulation 
requires adequate resolution of not only the 
support geometries, but also any 
anticipated induced flow features such as 
juncture vortices, or even localized flow 
separation. At a minimum, the support 

structures will influence the slot exit flow Figure 14. NTF air station operational requirements of 
quantities and may well dominate the model blowing system for different slot heights, To=-50°F, 

MAC=17.4”, b=44”, b SL o T /b=97.7%, Po=5 atm, M=0.2. 


i 5.00 
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internal flow field. Thus, documentation of the slot exit properties across the entire span of the wind tunnel 
becomes crucial. 

The importance of documenting the approaching upper surface boundary layer was also discussed in 
the 2-D CFD section, but will be expanded upon in the context of adding cross flow, due to either wing 
sweep or leading edge flow separations. Depending on the location and orientation of the blowing slot, it is 
possible that a highly three-dimensional mixing and entrainment process will occur just downstream of the 
blowing slot. In some applications, this may be the desired effect, but will only complicate off-body flow- 
field measurements necessary for proper CFD code validation. Also, the addition of flap side-edge and 
wing-tip vortices further complicate the flow structure near the model. Obtaining accurate experimental 
measurements for these types of flows will be a taxing exercise for state-of-the-art measurement 
techniques, and it will likely require several techniques to be used in concert to achieve quantifiable 
uncertainties. 

The following CFD examples highlight the unstructured flow solver USM3D being used to examine 
the effects of blowing at transonic and low-speed conditions. Representative Mach number contours shown 
in Figure 15 are for the modified NASA TMA-0712 airfoil. The Mach number is 0.78, the lift coefficient 
is 0.70, and NPR=1.25. The upper surface normal shockwave is just upstream of the blowing slot exit, and 
the blowing has been found to be effective in reducing the drag by reducing separation near the trailing 
edge of the airfoil. It should be noted that for the 3-D semispan model, this would correspond to a free 
stream Mach number of approximately 0.85. Additional studies are underway to optimize the low-speed 
high-lift components, depicted in Figure 16, for a lift coefficient of approximately 6.0. During the first 
phase of testing, the semi-span model will be fitted with a conventional leading edge slat. This will be 
replaced in subsequent testing with the leading edge blowing slot. 




Figure 15. Mach contours for the NASA TMA- 
0712 airfoil, M=0.78, NPR=1.25. 


Figure 16. Mach contours for the high-lift 
configuration for the NASA TMA-0712 airfoil, 
M=0.20. 


- 11 - 


Subsystem Case 

A Hybrid Wing Body (HWB) model with blown flaps (see Figure 17) was tested in the NASA Langley 14- 
by-22-Foot Subsonic Tunnel (14x22) at low speeds (q = 30 psf, = 0.143). The high-lift investigation 
was conducted as part of a 3 -way collaboration between NASA LaRC, the Air Force Research Laboratory 
(AFRL), and the Northrop Grumman Systems Corporation (NGC) to advance state-of-the-art technologies 
for CESTOL aircraft. The enabling tools and technologies for CESTOL aircraft are dual use in that they 
are applicable for both military and civil missions. From the perspective of NASA’s Subsonic Fixed Wing 
Project, the wind tunnel testing of the HWB model has provided unique databases to validate 3-D CFD 
tools for CESTOL efforts, particularly relating to the development of high-lift wings with advanced CC 
technologies. 



Figure 17. HWB model in the NASA Langley 14-by-22-Foot Subsonic Tunnel. 

The initial HWB test in the 14x22 used pressure sensitive paint (PSP) to complement pressure 
measurements from surface orifices. In addition, off-body flow field measurements were acquired with a 
7-hole-probe rake survey system. Although only a reduced number of nozzle pressure ratios and 
geometrical settings (i.e., angles of attack of 0° and 10°, flap deflections of 0° and 60°, with and without the 
slat) were possible in the first tunnel entry, the experiment revealed several “lessons learned” issues and 
challenges associated with the PSP and 7-hole-probe measurements that need to be addressed for future 
tests. 

PSP is an advanced pressure measurement technique that one would naturally consider for wind-tunnel 
testing of a full-span model because it can cover a substantial amount of surface area. PSP has 
demonstrated that it can provide reasonable visualization of surface pressure distributions globally. 
Nevertheless, the temperature gradients across the model due to a slow heat transfer associated with the jet 
plenums are a major challenge to overcome for the PSP technique. The plenums were preheated to the 
operational flow conditions prior to acquiring PSP data. While this reduced the temperature effect, the 
process was time consuming. One solution being developed for future testing is the use of a biluminophore 
PSP formulation to minimize the effect of temperature gradients. 

The off-body measurements were obtained using standard 7-hole probe techniques. Eight 7-hole probes 
were mounted on the rake head connected to a probe positioning traverse system that was used to survey 
the wake of the left wing. All U, V, and W velocity components are in the tunnel coordinate system instead 
of model coordinate system. Figure 18 shows typical 7-hole probe results capturing complex and 
interesting flow features in the cross plane between the inboard flap edge and the outboard wing’s trailing 
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edge. The two examples are the no blowing and blowing cases (NPR of 1.0 and 1.37, respectively) for the 
model with the slat installed, a flap deflection of 60°, and an AOA of 10°. The data are presented in terms 
of streamwise mean velocity (U) contours and velocity vectors in the Y-Z plane. 


Run 2229, NPR = 1.0 wife 
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Figure 18. U-velocity contours and velocity vectors in the Y-Z plane 
(Moo = 0.143, 60° flap and slat on, AOA = 10°). 

Two swirling regions downstream of the outboard flap edge and the inboard/outboard flap juncture are 
highlighted with velocity vectors for the blowing case (NPR = 1.37) shown in Figure 18. As expected, the 
increased lift associated with the 60° blown flap significantly enhanced the downwash velocity and the 
swirling. The stronger swirling is associated with a flap-edge-induced outboard vortex that has 
significantly more velocity deficit in its core than the weaker one downstream of the flap juncture. When 
viewing toward the upstream direction, the swirling motion is clockwise for both vortices. The outer edge 
of a third vortex with counter-clockwise rotation is caused by the inboard edge of the inboard flap and can 
be observed near the right edge of the plots. 

White blanks in the middle of a vortex or wake for all cases indicate bad data (erroneous probe total 
pressure) due to high angularity in the local flows. Flows near the core of a wake vortex or near the off- 
body re-circulating regions often had flow angles greater than 75° with respect to the probe axis, which are 
outside the measuring capability of the 7-hole probes. 

The 7-hole probe technique has demonstrated its capability of identifying complex and interesting 
flow features in the wake downstream, such as velocity-deficit and swirling-flow regions, for CFD code 
validation and for enhancing the understanding of the complex flow physics. In addition to the 7-hole 
probe measurements, a large field of view Particle Image Velocimetry (PIV) was used during the second 
wind tunnel entry in the wake of the model. The PIV system had some advantages such as non-intrusive 
and simultaneous acquisition of a plane of data. However, PIV data does not provide pressure information 
that can be correlated directly to the 7-hole probe data. It is also recognized that the particle density in the 
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strong vortex cores was too low to accurately resolve the velocities and that high tunnel temperature also 
influenced the particle density. Using both 7-hole probe and PIV data with uncertainty estimations is 
useful for CFD validation. 

Documentation of the slot flow and boundary layer using a high resolution hot-wire survey is the best 
way to provide the blowing boundary conditions for CFD. In short, several different measurement 
techniques (i.e., force and moment, pressure orifices, PSP, 7-hole probes, PIV, hot-wires, Preston tubes, 
etc.) must be used in concert to provide a full database for CFD validation. 

Complete System Case 

The advances realized through this CFD validation effort can be applied to integrated CC technologies 
applied to a CESTOL aircraft. Application of CC to a HWB concept, such as shown in Figure 19, will 
require the integration of several advanced systems. The development of a large-scale wind tunnel test or 
flight test for such an aircraft will require collaboration with industry and other government agencies. 



Figure 19. Sketch of a hybrid wing body configuration with 
circulation control technology applied to the outboard wing 
sections. 


Concluding Remarks 

This paper highlights on-going research efforts at NASA concerning circulation control that emphasize 
CFD validation and the creation of benchmark experimental databases to advance our knowledge base and 
prediction capabilities. The systematic CFD validation approach that is described in this report is 
consistent with NASA’s foundational flow physics research for near-term conventional aircraft applications 
and longer-term hybrid wing-body configurations. The measurement and validation challenges for small- 
scale circulation control experiments are being pursued through the approach described in this paper. The 
circulation control examples shown in this paper highlight recent advances in CFD and experimental 
methods. These examples set the stage for developing wind tunnel databases related to supercirculation 
that will be required for CFD validation. The planned high-Reynolds-number experiments in NTF will 
provide a unique data set that will enable researchers to separate out the effects of the different blowing 
parameters so that they can better understand the complex physics of advanced circulation control 
applications. The outcomes of these efforts are expected to improve both experimental and CFD 
capabilities related to powered-lift and circulation control. 
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